Spin-1/2 XX Z diamond chain within the Jordan- Wigner fermionization approach 
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The spin-1/2 XX Z diamond chain is considered within the Jordan- Wigner fermionization. The 
fermionized Hamihonian contains the interacting terms which are treated within the Hartree-Fock 
approximation. We obtain the ground-state magnetization curve of the model for some particular 
cases and compare the results with the exact diagonalization data for finite chains of 30 spins and 
known exact results. We also analyze the vahdity of the suggested approximation. 
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The spin-1/2 XX Z diamond chain as a quantum frus- 
trated system is an interesting topic for the theoretical 
research, since it exhibits many phenomena related to 
interplay of quantum fluctuations and competing inter- 
actions. The exact results for this model are available 
only for limited cases [l|, Q , and the earlier studies used 
mainly the numerical methods 0411 • We suggest an ana- 
lytical method based on the fermionization of the initial 
spin-1/2 model. The present work is an extension of our 
previous study of the XX diamond chain using the 
Jordan- Wigner transformation and the Hartree-Fock ap- 
proximation. The goal of the paper is to explain the 
properties of the frustrated model on the fermionic lan- 
guage for a more general XXZ model. We accompany 
our approximate calculations by the results of the exact 
diagonalization and also compare them with the partic- 
ular cases where the exact results are available [l|, . 

We consider the quantum spin-1/2 XXZ model on the 
generally distorted diamond chain (sec Fig.l in Ref@) 
with the following Hamiltonian: 
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(1) 



Here, s^j (a = x,y,z) are the Pauli spin-1/2 operators 
with the first index corresponding to a sublatticc and the 
second index to a cell, Jp = J| = Jp > 0, = AJp > 0, 
A is the interaction anisotropy, and h is the external 
magnetic field (we set g^iB = !)• Using the spin raising 
and lowering operators s^^ = s^i ztiSp ^ one can rewrite 
the xy-part of the Hamiltonian as a quadratic form of 
the mentioned operators. The zz interaction, due to the 



relation sf 



1/2, leads to the product of four 



Following the procedure described in Ref|6|, we specify 
the Jordan- Wigner transformation in the form: ci^i = 

11^=1 ntl-P-J.iSr,;' C2,; = Ul=lU^i=[Pq,iPl,lSl,l, C2,l = 

nl=iUt\P,-^Pi,iP2.isli, where P,, = (-25^,) = 
exp{iTTc^ iCq,i) is the Jordan- Wigner factor. Here new 
operators Cp,i satisfy the Fermi commutation relations. 
The fermionic expression for the z-componcnt of the spin 
operator can be easily obtained as s; 



q,l 



1/2. 

Therefore, the zz-part of the Hamiltonian contains the 
four-fermion term. The xy-part of the Hamiltonian in 
terms of new Fermi operators is as follows: 
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+ J3(c+,P2,;C3,/ + C+;P3,;Ci,; + i) + h.C.) 
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spin-lowering and raising operators. 



Note, that the transformed Hamiltonian contains the 
fcrmion interaction for terms proportional to the J3 cou- 
pling which is usually assumed to be smaller than J2. We 
should also note the relation between the spins- 1/2 and 
spinless fcrmions: the spin-down (-up) state corresponds 
to the empty (filled) fermionic state; the action of the 
fcrmion creation and annihilation operators is analogous 
to the action of spin raising and lowering operators in spin 
language. The resulting Hamiltonian represents the in- 
teracting Fermi gas. To proceed we use the Hartree-Fock 
approximation where all interacting terms are factorized 
preserving all pair correlations between nearest neigh- 
bors of type {CpiCq^m}- Thus, the Hamiltonian becomes 
a quadratic form in terms of Fermi operators. It can be 
diagonalizcd using Fourier and Bogolyubov transforma- 
tion, and the thermodynamics and fermionic correlation 
functions are easily found. However, it depends para- 
metrically on the unknown contractions {cyiCq,m) which 
have to be found self-consistently. 
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FIG. 1: The ground state magnetization versus external field 
for a) Ji = J3 = 1, J2 = 1.75; b) Ji = 1, J2 = 1.25, J3 = 0.45. 
Step-like dashed lines represent the exact diagonalization data 
for 30 spins, solid lines correspond to the approximated re- 
sults, dotted lines indicate the magnetization jumps. 

Solving the self-consistent equation we found that the 
dimer- monomer ground state for the symmetric dia- 
mond chain [l[ is recovered. Indeed, if J 2 > 2Ji and 
the external field < h < JiA -I- J2(l + A)/2, we 
obtain the following solutions for the elementary con- 
tractions: {c'ljCij) = 1, (cJ,C2,;) = (c^jC^^i) = 1/2, 

{ct.iC2,i) = {ct^iC3,i) = 0, (cJ;C3,/) = -1/2. The ground 
state of the corresponding fermion model is as follows 
\GS) = Uict/c+i - c+,)/a/2|0) where |0) denotes the 
empty state. Using the relation between fermionic and 
spin states one can see that it corresponds to the dimer- 
monomer state in spin language (see also Ref0). The 
results for the symmetric diamond chain below the dimer- 
monomer limit is shown in Figllji. We observe an excel- 
lent agreement with the exact diagonalization data ex- 
cept for small fields where our mean- field- type approxi- 
mation produces the non-zero magnetization. As an ex- 
ample of the distorted chain, we have chosen the param- 



eter set used previously for azurite In FiglT|D we ob- 
tain a good agreement with the XX limit. Particularly, 
the magnetization curve shows the 1/3-plateau and zero- 
magnetization in zero field. However, for the isotropic 
Heisenberg interaction the discrepancy between the ap- 
proximate and exact results raises up quickly below the 
upper critical field. The exact diagonalization data for 
A = 1 seem to show a cusp in the magnetization curve 
at 2/3 of the saturation magnetization. It might be the 
sign that the magnetic cell is doubled for this model. As 
it was discussed in 0, to describe this kind of behavior, 
it is necessary to consider also non-uniform elementary 
contractions. Another drawback of the Hartree-Fock ap- 
proximation is the artificial jumps of the magnetization. 

We have also examined our method by comparison 
with the exact results for the Ising-Heisenberg diamond 
chain obtained by means of decoration-iteration proce- 
dure ■ It is the special case of the anisotropic diamond 
chain where the spins on the vertical bond are coupled 
by the Heisenberg interaction whereas all other couplings 
arc of Ising type. The model ([T]) corresponds to the model 
considered in if we put J^'^ = J3 = 0, J2 = AJ2, 
Ji ~ Ji = Ji, Ji ^ J^- Depending on the ratio be- 
tween the interaction couplings, anisotropy and external 
field the system may stay in different phases: |FRI) = 
ri; I ;i,/>l t2,/t3,/>, |FRU) = Hi I ti.iXI t2,/;3,/>-| ;2,;t3,/ 
))/V2, |SPP) = Hi I UiUi). Within our Hartree- 
Fock approach it is possible to recover the ground state 
phase diagram and all phase boundaries. The men- 
tioned states in the fermion representation have the form: 
|FM) = c+;<ilO), |FRU) - Ui <,(4z " <,)/V2|0), 
|SPP)=nz<;<;C+|0). 

To conclude, we have considered the approach based 
on the Jordan- Wigner fermionization and subsequent 
Hartree-Fock approximation for the spin-1/2 anisotropic 
diamond chain. We have revealed that such an approach 
recovers the exact results for the phases characterized 
by short-range correlations as, for instance, the dimer- 
monomcr phase. It also provides a good description of 
the ground state magnetization for the symmetric dia- 
mond chain. However, for the distorted diamond chain 
the inclusion of the zz interactions between spins may 
lead to the qualitative change of the magnetization curve 
which cannot be explained within the current approxi- 
mate method. 
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